Direct imaging of topological edges states with cold atoms 
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Detecting topological order in cold-atom experiments is an outstanding challenge, the resolution 
of which offers novel perspectives on topological matter. In material systems, unambiguous signa- 
tures of topological order exist for topological insulators and quantum Hall devices. In quantum Hall 
systems, the quantized conductivity and the associated robust propagating edge modes - guaranteed 
by the existence of non-trivial topological invariants - have been observed through transport and 
spectroscopy measurements. Here, we show that optical-lattice-based experiments can be tailored 
to directly visualize the propagation of topological edge modes. Our method is rooted in the unique 
capability for initially shaping the atomic gas, and imaging its time-evolution after suddenly remov- 
ing the shaping potentials. Our scheme, applicable to an assortment of atomic topological phases, 
provides the first method for imaging the dynamics of topological edge modes, directly revealing 
their angular velocity and spin structure. 



The integer quantum Hall (QH) effect revolutionized 
our understanding of quantum matter, revealing the ex- 
istence of exotic phases not described by the standard 
theory of phase transitions [1]. In this phenomenon, the 
Hall conductivity is quantized, anh/e 2 = v E Z (where e 
is the electron charge and h is Planck's constant), when- 
ever the Fermi energy resides in an energy gap. The 
integers v are related to topological invariants - Chern 
numbers - that are associated with the bulk energy bands 
[1, 2]. Their topological origin guarantees that the Chern 
numbers are constant as long as the bulk gaps remain 
open, explaining the signature plateaus in the Hall resis- 
tivity, present when external parameters, such as mag- 
netic fields, are varied. Moreover, a holographic princi- 
ple stipulates that a topologically-ordered bulk gap, with 
topological invariant z/, necessarily hosts v propagating 
modes localized on the samples edge [3] . These topologi- 
cal edge states are chiral - their motion has a well-defined 
orientation - inhibiting scattering processes in the pres- 
ence of disorder. 

In condensed matter physics, direct observations of 
edge states remain relatively rare. A first signature was 
obtained from magnetoplasmons created by pulsed volt- 
ages [4]. Another evidence arises from edge transport 
in engineered Aharonov-Bohm interferometers with QH 
systems [5]. By contrast the "routinely used" spectro- 
scopic reconstruction of mid-gap states [1] is consistent 
with the expected topological band structure, but does 
not prove their chiral nature. 

Cold atoms trapped in optical lattices and subjected 
to synthetic gauge fields [6, 7] are an ideal platform for 
realizing topological insulating phases. Making topology 
manifest in experiments, however, is a fundamental chal- 
lenge. In this context, transport experiments are conceiv- 
able [8] but technically demanding. Existing proposals 
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FIG. 1: Experimental scheme and general strategy. 

a, Trapped cold atomic fermions move on a square optical 
lattice in the presence of a synthetic uniform magnetic flux 
<I>. Two repulsive potentials, initially piercing holes in the 
atomic cloud, are suddenly removed at time t = 0. At all 
times, atoms are confined by an additional circular poten- 
tial. We assume in this paper that the confining barriers are 
perfectly sharp, and discuss the case of smoother potentials 
in the Appendix E. b, The system is initially prepared in a 
quantum Hall phase: chiral edge states propagate along the 
edges determined by the repulsive walls and the external con- 
finement. After releasing the walls, the edge states tend to 
propagate along the Fermi radius determined by the circular 
confinement: they encircle the initially vacant regions. 



for measuring topological invariants [9-13] have experi- 
mental drawbacks and can only be applied to very spe- 
cific configurations. Likewise, detecting topological edge 
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states [14-19], e.g., through light scattering methods [20], 
require complicated manipulations to separate the small 
edge- state signal from the bulk background [17]. 

Here, we introduce a simple method to directly 
visualize the propagation of topological edge modes, by 
studying the time evolution of an atomic QH system 
after suddenly releasing constraining walls, see Fig. 
1. We show that the movement of the chiral edge 
states, encircling the initially vacant regions, is then 
directly visible in the atomic density. This allows us 
to evaluate the edge states angular velocity and to 
provide an unambiguous signature for the existence 
of chiral edge modes in the atomic system (Fig. 1). 
Our method is straightforward and insensitive to 
experimental imperfections: it only relies on a large 
initial occupancy of edge states. Crucially, our method 
requires that the edge states contribution to the density 
remains spatially separated from the bulk, which can be 
realized by populating a dispersionless bulk band with 
non-zero Chern number. We present several detection 
techniques, applicable to both flat and dispersive bands, 
that demonstrates the universal applicability to atomic 
systems with propagating edge modes. 

We consider a two-dimensional optical lattice filled 
with non-interacting fermions, subjected to a uniform 
synthetic magnetic flux $ [21, 22], and confined by a 
circular potential, V con f(r) = Vo(r/ro) 7 . In experiment, 
V con f(r) can be made nearly arbitrarily sharp (7 — >• 00) 
[23, 24]; this configuration is of particular interest for 
our scheme, as demonstrated in the Appendix E. The 
resulting system realizes the Hofstadter model [25] with 
second-quantized Hamiltonian 

(i) 

m,n 

cj^ n describes the creation of a fermion at lattice site 
x/a = (m, n) G Z 2 ; J is the tunneling amplitude; and 
we take the lattice period a as our unit of length. This 
model has a topological band structure [1, 2]: When 
= p/q G Q, the bulk energy spectrum splits into q 
subbands [25], each associated with a non-zero Chern 
number [2] . This guarantees the existence of robust edge 
states in the bulk energy gaps [3]. These edge states 
are chiral in the sense that they propagate along the 
Fermi radius Rp (i.e., the edge of the atomic cloud de- 
limited by the confining potential V con {) with a definite 
orientation of propagation. It is convenient to represent 
such non-trivial spectra by diagonalizing the Hamilto- 
nian (1) on a cylindrical geometry [3], see Fig. 2. This 
picture shows a clear separation of the bulk and edge 
states dispersions that survives in the experimental cir- 
cular geometry produced by V con { (Appendix E and Refs. 
[17, 26]). In the following, we specifically study the con- 
figurations $ = 1/3 and <£> = 1/5 and set the Fermi en- 
ergy = —1.5 J inside the lowest bulk energy gap, see 
Fig. 2. In both these situations, the lowest energy band is 



associated with the Chern number v = — 1, which guar- 
antees the occupancy of a single edge mode with same 
chirality sign(^) = (— ). These two configurations differ 
significantly in that the occupied bulk band is nearly flat 
(dispersionless) in the case <I> = 1/5, while it is dispersive 
for $ = 1/3, see Fig. 2. 
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FIG. 2: Bulk and edge states spectrum: dispersive 
vs flat bands. Energy spectrum E(k y ) as a function of the 
quasi- momentum k y for a, <1> = 1/3 and b, $ = 1/5, obtained 
by diagonalizing the Hamiltonian (1) on a finite cylinder di- 
rected along the x direction, with Konf = 0. The projected 
bulk bands E(k x ,k y ) —> E(k y ), shown in blue, are separated 
by large gaps of order ~ J. The red dispersion branches that 
are visible within the bulk gaps correspond to propagating 
modes that are localized on the opposite edges of the cylin- 
der. When the Fermi energy is set within the first bulk gap, 
a single edge mode is populated on each edge of the cylin- 
der (the lowest bulk band corresponds to the Chern number 
v — — 1 for $ = 1/q). When considering the circular geometry 
realized in an experiment (V CO ni / 0) and setting E? = —1.5 J, 
one is guaranteed that a single edge mode will be populated 
since the Chern number v does not rely on the specific ge- 
ometry used [17, 26]. Note that in the case $ = 1/5, the 
lowest energy band is characterized by the tiny flatness ra- 
tio, / = W/A w 0.04, where W(A) denotes the first band 
(gap) width; in this topological quasi-flat band configuration, 
the populated edge states are expected to propagate more 
rapidly than the bulk. 

Our scheme (a) demonstrates the existence of prop- 
agating modes that are localized close to the Fermi 
radius, and (b) identifies their angular velocity 0. 
To achieve this goal, we consider a geometry which 
constrains the QH system within two regions of the 
trap, as sketched in Figure 1, resembling a bat in flight. 
This initial "bat" geometry is shaped by a pair of 
sharp potential walls Vhoie = Vhoie 1 + Vhoie 2 defined 
by (x ± r /2) 2 + (y/V2) 2 < Tq/4, creating holes in the 
density distribution (see also Appendix F). In the bat 
geometry, we set the Fermi energy within the lowest 
bulk gap E F = —1.5 J and suddenly remove Vhoie a t time 
t = 0. We then study the dynamics of the atomic density 
with all other parameters unchanged. The bat shape is 
optimized for visualizing the time-evolving chiral edge 
states in the density p(x, t) for t > 0, see Fig. lb. In the 
following, we discuss how this "wall-removal" strategy 
can be exploited to reveal the edge states properties, as 
they progressively encircle the initially empty regions in 
a chiral manner. 

Figure 3a shows the time-evolving density distribution 
p{x,t) for $ = 1/3. This example highlights the im- 
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portance of the bulk band structure, as it demonstrates 
the drawbacks encountered when preparing a system 
with dispersive bulk bands, see Fig. 2a. The time 
evolution in Fig. 3a, illustrates two main effects: (1) 
the progressive encircling of the holes by particles at the 
system's radial-edge (with localization length ~ a), and 
(2) the rapid filling of the holes by bulk states, see also 
Fig. lb. Once t « 10 — 20ft/ J, the atomic cloud's initial 
bat shape has become cyclonic, directly demonstrating 
the edge states' chirality. To provide further insight, 
we separately calculated the contribution stemming 
from the initially populated edge states, p e d ge (x,t), see 
Appendix A. In the corresponding Fig. 3b, we observe 
that the edge states, which propagated along the edges 
delimited by the bat potential at t < 0, become localized 
along the circular edge at r = Rp-, and that they follow 
a "chiral" motion. These edge states remain localized on 
the edge for very long times, and only slightly disperse 
into the bulk of the system, as can be anticipated 
from the small wavefunction overlap between edge and 
bulk states. Figure 3a emphasizes the problematic 
(non-chiral) filling of the holes by the many dispersive 
bulk states. The speed at which this filling occurs is 
to be compared with the circular motion of the edge 
states, which can be estimated from the group velocities 
v = (l/h)dE/dk associated with the bulk and edge 
states, see Fig. 2a. Additional interference takes place 
within the system, leading to small but visible patterns 
in the density. Figure 3b, shows that these patterns 
mainly originate from interferences between bulk states. 

For $ = 1/5 and E F = —1.5 J, the dispersionless (flat) 
bulk band represented in Fig. 2b is totally filled, and 
most of the edge states lying in the first bulk gap are 
populated. The corresponding time-evolving density 
p(x, £), depicted in Fig. 4a, is radically different than for 
<I> = 1/3, cf. Fig. 3a. For <l> = 1/5, the edge states encir- 
cle the initially forbidden regions, largely unperturbed 
by the now motionless bulk. The dispersionless nature 
of the bulk states is further illustrated in Fig. 4c, which 
shows the evolution of pbu\k(&,t) — p(x,t) — p e & ze (x,t). 
The initial bulk states are immobile for times ~ 30ft/ J. 
In Fig. 4b, the evolution of the initially populated 
edge states p e dge(x,t) shows an interesting behavior: 
edge states with energies close to the flat bulk band 
are almost dispersionless, and remain localized on the 
edges delimited by the bat potential. Besides, the many 
dispersive edge states at higher energies encircle the 
holes in a clear and chiral manner. The instructive 
dynamics of the density p(x,t), which is due to the 
clear separation of the edge and bulk states during 
the evolution, can also be understood by studying the 
population of the single-particle eigenstates (Appendix 
A). Moreover, our method is highly robust against 
perturbations in the density (or equivalently in the 
Fermi energy, ~ —1.5 J + 5), as it only relies on the 
occupation of dispersionless bulk states and sufficiently 
many edge states. We verified that the edge state signal 
is unambiguous when a high energy (dispersive) band 
is initially filled (Appendix D). We also found that this 
signal is immune to an additional random potential 
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FIG. 3: Evolution of the spatial densities after re- 
leasing the walls, a, the spatial density p(sc,£), and b, the 
contribution of the initially populated edge states p e d ge (#,£)• 
In all the figures, $ = 1/3, Ef = —1.5 J and we considered in- 
finitely sharp circular confinement (ro = 27a) and ellipsoidal 
walls Vhoie- The total number of particles is iVp ar t = 210, 
while the number of initially populated edge states is iV e dge ~ 
80. 
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FIG. 4: The topological quasi-flat band configura- 
tion, a, the spatial density p(x,t), b, the contribution of the 
initially populated edge states p e dge(x,t), and c, the contri- 
bution of the initially populated bulk states pbuik (#,£)• In 
all the figures, $ = 1/5 and Ef — —1.5 J. The total num- 
ber of particles is iVp ar t = 146, while the number of initially 
populated edge states is N e d ge ~ 64. Note the dispersionless 
nature of the occupied bulk states, which highly improves the 
detection of the edge state signal. 



|Kand(#)| i$ J mimicking disorder, in agreement with 
the chiral nature of the propagating states (Appendix 
G). Thanks to the topological quasi-flat band config- 
uration, the edge states propagation can be visualized 
on long time scales, without being affected by the bulk 
dispersion. For = 1/5 and a typical system size 
Rp ~ 100a, we find ~ O.OU/ft, which would require 
a realistic time t ~ 70ft/ J to observe the 7r/4-rotation 
undergone by the edge states. 



4 



a £p = 


= p(x,t; $ 


= +1/3) - 


- p(sc,£; $ = 


-1/3) 






% 














\ 


J 

21H/J 






A9H/J 






70H/J 


^ Pedg 


s (x,t;$ = H 


-1/3) 












1 
















V 


21h/J 






49H/J 






/ 

70H/J 




} (a?,t;$ = - 


-1/3) 














21/i/J 






4Qh/J 






70H/J 



0.2 



-0.2 
0.2 



I 





0.2 



FIG. 5: The opposite-flux method for dispersive 
systems, a, Evolution of the difference 5p = p(x,t;<fr = 
+1/3) — p(x,t;<f> = —1/3), for the same configuration as in 
Figure 3. This method yields a clear manifestation of the edge 
states, 5p « p e d g e(+^ ) ) — p e dge(— 3>), by eliminating the unde- 
sired contribution of the many bulk states, b, The edge- states 
contribution p e dge(x,t) for $ = +1/3, and c, for $ = —1/3. 
In the central column, we note the vanishing of the signal 
Sp(x) ~ that occurs at time t* ~ 49ft/ J, indicating that the 
edge states angular velocity is ~ 0.03 J/ h for Rf = 27a and 
7 = oo (also see Appendix B). 



We now describe two methods for isolating the edge- 
states contribution p e dge from that of the bulk states, 
useful for systems with dispersive bulk bands. The first 
method consists in performing two successive experi- 
ments, using the same apparatus and parameters, but 
with opposite fluxes and — The difference between 
the two images 5p(x, t) = p(x, t; +<£) — p(x, t\ — $), lacks 
the non-chiral contribution of the bulk states and is sim- 
ply given by Sp w p e d g e(+$) - Ped g e(-$), see Fig. 5. 
For our bat geometry, p e dge(+^) ~ Pedge( — when 
the edge states have undergone a rotation of 6 = tt/2. 
This determines the time t* = it / (26) when the signal 
5p(x,t*) disappears, giving the angular velocity of the 
edge states (Appendix B). This situation is illustrated in 
Fig. 5 for the $ = 1/3 "dispersive" case, where we find 
5p(x,t* ~ 49ft/ J) « 0, in good agreement with the an- 
gular velocity e « 0.03J/ft of the populated edge states 
[17]. We verified that slight differences in the filling [e.g., 
£ F ($+ = +1/3) - £ F ($_ = -1/3) ± 0.1J], for example 



due to finite temperature effects between the two succes- 
sive experiments, or variations in the flux (e.g., 3> + = 1/3 
and <£_ « —$+±0.01), do not significantly affect the sig- 
nal 5p(x,t), highlighting the robustness of this method 
against experimental imperfections. 

The second method aims to efficiently reduce the bulk 
dispersion by suddenly lowering the potential walls Vhole 
at t = 0, instead of removing them completely This 
operation can be achieved in such a way that only the 
edge states with sufficiently high energies are allowed 
to propagate, while leaving the bulk states away from 
the holes. This "edge-filter method" can be realized 
by setting the Fermi energy within the first bulk gap, 
and then suddenly lowering the potential Vhole to the 
value Vjfde ~ w at * = °> where W is the width of the 
lowest bulk band. The great efficiency of this method is 
presented in the Appendix C for the case <1> = 1/3. 



Our results are stable when the spacious ellipsoidal 
walls are replaced by small perturbative potentials 
(Appendix F). Moreover, using a gaussian potential 
instead of the infinitely sharp potential Vhole gives 
equivalent results (Appendix E). In contrast, a sharp 
external confinement V con f(r) ~ r 1 with 7 ~ 00, 
is desirable for detecting the edge states motion in 
reasonable times, since the edge-state velocity dramat- 
ically decreases as the potential is smoothened [17]. 
Furthermore, the radii R e at which the edge states (j) e 
localize are energy-dependent R e ~ yfel, leading to a 
spatial broadening of the edge state signal for finite 7 
(Appendix E). The schemes introduced in this work to 
reduce, or even eliminate, the irrelevant contribution of 
dispersive bulk states can be applied to a wide family 
of topological atomic systems, such as the promising 
Haldane-like optical lattice [9, 11] and fractional QH 
atomic gases [27, 28]. Finally, our method can be directly 
extended to visualize the propagation of Z 2 topological 
(spin-polarized) edge states, both in 2D [15, 29] and 3D 
[30], by using standard spin-dependent imaging methods. 
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Appendix A: The time evolution and observable 
quantities 

The system is prepared in the ground-state of the ini- 
tial Hamiltonian, 

(Al) 

where the potential Vhoie describes the walls initially 
present in the trap (see main text), and where H is 
given by Eq. (1). We denote the number of avail- 
able sites before and after removing the walls by 
ng ites and n S it es , respectively, and we define the ratio 
(1 - rj) = n° tes /n sit es. When V ho}e = 0, the total 
number of sites within the trap is approximative ly 
given by the area n s i te ~ 7r(ro/a) 2 , while the outer 
circular edge contains about n e dge ~ 2irro/a lattice sites. 
The holes in the density created by V^oie correspond 
to n ho ie = ^sites - ^sites ~ 7r(r /a) 2 / V2 vacant sites, 
leading to the large ratio r] = nhoie/^sites ~ l/\/2- 
For ro = 27a, the system initially contains ng ite ~ 700 



sites and the number of sites delimiting the edge of 
the bat is n® dge « 320 sites. Thus, the spacious holes 
used in our calculations lead to a large edge/bulk 
ratio. For $ w p/q £ Q, and initially setting the Fermi 
energy in the lowest bulk gap, leads to the filling factor 
v° = ^part/^sites ~ After removing the walls 

Vh i e , the filling factor is reduced to the smaller value 

V = iVpart/^sites - (1 - Tj) / q <C zA 

The groundstate of Hamiltonian (Al) is written as 

i*o>= n /«i >> ( A2 ) 

E a <E F 

where the operator /t creates a particle in the single- 
particle state \xa), with energy E a located below the 
Fermi energy Ep. Here {\x a ),E a } represents the com- 
plete set of single-particle eigenstates and eigenvalues sat- 
isfying the stationary Schrodinger equation 

H \Xa) =E a \ Xa ). (A3) 

We are interested in the time evolution of the spatial 
density p(x,t) after removing the walls Vhoie at t = 0. 
The evolution of the single-particle states \x a ) is then 
entirely governed by the Hamiltonian H. It is therefore 
convenient to introduce the eigenstates and eigenvalues 
{\4>\)-> e \\ corresponding to the Hamiltonian H, 

H\<f> x ) = e A |0 A >. (A4) 

We then introduce the notation \ X a(t)), which corre- 
sponds to the time evolution of the initial state \ Xa ) ac- 
cording to the Hamiltonian H, 

\x a (t)) = J2(^\x°) e ~ itxt/H \^)- ( A5 ) 

A 

The spatial density p(x, t) at time t is given by 

P (x,t) = (*(t)|4,nCm,n|*W) 

= <*0|4,n(*)2m,nWI*0), (A6) 

where |^o) is the initial ground state in Eq. (A2). Here, 
we introduced the Heisenberg representation of the field 
operators 

A(t) (+) = p iHt/h*ti) p -iHt/h 

= £ x M(*,i)/(t), (A7) 

a 

where Xa(x,t) = (x\ X a(t)) is given by Eq. (A5). From 
Eqs. (A2),(A6),(A7), one finds the simple result 

p(x,t)= \Xa(x,t)\ 2 , (A8) 

E a <E F 

namely, that the particle density p(x,t) is entirely gov- 
erned by the time-evolution of the initially occupied 
single-particle states. The time evolution of the atomic 
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cloud, after releasing the walls Vhole at £ = 0, can there- 
fore be numerically evaluated through a direct diagonal- 
ization of the Hamiltonians H and Hq. In our study, 
a crucial aspect consists in identifying the regimes for 
which the edge states propagating around the initially 
forbidden regions provide a clear signal, which is not per- 
turbed by the many bulk states. It is therefore desirable 
to separately evaluate the contributions of the initially 
populated bulk and edge states. We introduce the corre- 
sponding quantities 

PedgeO,£) = IXeO,£)| 2 , 
E e <E F 

Pbulk(#, t) = p(x, t) - Pedge(x, t), 

where the sum ^2 E <Ef is restrained to the populated 
edge states with energies E e located within the bulk gap. 

In this study, the time evolution is chosen to be en- 
tirely dictated by the Hamiltonian H, and thus, it is 
non-dissipative: after releasing the walls, the total en- 
ergy of the system is constant and is given by 

So = <*o |#|*o) « Yl E «> ( A9 ) 

E CX <E F 

where we considered the approximation (^o\H\^o) « 
(^o |-ffo I ^o)> which is valid for Vhole sufficiently abrupt. 
For example, in the situation illustrated in Fig. 4, 
the energy released after removing the walls is about 
J/3. Therefore, in our calculations, the many-body state 
\^(t)} never reaches the ground state |^gs) of the final 
Hamiltonian H, with energy 

£gs = <^gs|£|^gs). (A10) 

This final ground state, which differs from the initially 
prepared ground state |^o) in Eq. (A2), is expressed as 

A/p ar t 

i* GS ) = n a !i0>. (ah) 

A=l 

where creates a particle in the single-particle state 
|0 A ), with energy e A < e A +i- Here, 7V part = Y, Ea <E F 
is the total number of particles in the system, which is 
supposed to be constant (see Figs. 6a-b). In our non- 
dissipative framework, the probability V\(t) of finding a 
particle in the eigenstate \<j)\) is constant and inferior to 
one, as it is simply given by 

V x (t) = (*(t)\a{a x Mt)) = Yl l«>A|Xa}| 2 =cst<l, 

E a <E F 

^ (*cs|aiaA|*Gs) = 1 if A < iV part (otherwise 0). 

(A12) 

The populations T\ are illustrated as a histogram in Fig. 
6c, for the case $ = 1/5 and E F = —1.5 J. We find that 
when the Fermi energy is initially set within the first bulk 
gap, the population of high-energy dispersive bulk bands 
is highly limited during the evolution, which guarantees 
a clear spatial separation of the bulk and edge signal in 



this "topological quasi-flat band" configuration. The effi- 
ciency with which the initial edge modes Xe project to the 
final edge modes <j) e is further shown in Fig. 7, indicating 
the success of our scheme. For the situation illustrated 
in Figs. 6-7, we find that the number of populated edge 

States V edge = EAGedge states ~ 30 ~ ^edge/2, where 

A^edge ~ 64 is the number of populated edge states be- 
fore removing the walls (see main text). This result is in 
agreement with the fact that the edge delimited by the 
bat contains n® dge « 2xn e d ge sites, where n e d ge ~ 27rro/a 
is the number of sites delimited by the external confine- 
ment. 

In the main text, we discuss the non-dissipative evolu- 
tion of the atomic cloud after releasing the walls Vhole, ne- 
glecting thermalization processes. However, it is instruc- 
tive to estimate the energy loss that would be required 
to reach the ground state |^gs) 5 namely SE = Eq — Eqs- 
We propose to evaluate this energy difference for a con- 
figuration which is particularly relevant for our work (see 
main text), namely, a system exhibiting large initial holes 
and a flat (dispersionless) lowest bulk band. In the fol- 
lowing, we set the energy of the lowest bulk band equal 
to zero. The energy difference SE will generally be large 
in this "large hole/flat band" situation, since many bulk 
states become available in the lowest energy band after 
removing the walls, namely Eqs ~ [see Figs. 6a-b and 
below] . For this configuration, which is particularly suit- 
able for visualizing the edge states encircling the holes, 
the energy difference is thus approximatively given by 

6£* Yl E ^ ( A13 ) 

occupied edge states e 

where we considered the Fermi energy to be located inside 
the first bulk gap [see Fig. 6a-b]. For large holes r] = 
^hoie/^sites ^ 0, the final filling factor v <C 1/q such 
that the lowest bulk band will only be partially occupied 
when reaching the ground state I^gs)? and thus Eqs ~ 
0. From a rough geometric argument, we expect about 
A^edge ~ ttRf(2 + y/fj) / a(q — 1) available edge states in 
the lowest bulk gap of the initial system. Approximating 
the edge state branch as being linear inside the whole 
gap A, we find 

SE ~ -^-(R F /a)(l + 0.507), (A14) 

which corresponds to an energy difference SE ~ 
1.5(-Rp/a)J in the case $ = 1/5. In a typical exper- 
iment, Rf ~ 100a, which would require an important 
energy loss SE ~ 150 J to reach the ground state |^gs)- 
We verified that the estimated energy difference SE in 
Eq. (A 14) is in good agreement with a direct numerical 
evaluation of Eq. (A13), for the system illustrated in Fig. 
4 (i.e., $ = 1/5, r] « l/y/2 and E F = -1.5J). 

Appendix B: The opposite-flux method for 
dispersive bulk bands 

In this Section, we discuss how the difference 5p(x, t) = 
p{x,t\+§) — p{x,t\— <I>), obtained from two successive 
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FIG. 6: Ground states and single-particle states pop- 
ulation. Comparing a, the groundstate \*&o) of the initially 
prepared system (with holes) and b, the groundstate |^gs) 
of the final system (without holes). Filled and empty blue 
(red) dots represent the occupied and unoccupied bulk (edge) 
states, respectively, with occupation number = 1. Note that 

is constant. 



the total number of particles, iY pa rt 



*<£f' 

c The population V of the states \4>\), as a function of their 
energy e\, as established by the Fermi energy Ef = —1.5 J for 
the case illustrated in Fig. 4 (i.e., $ = 1/5, 7 = 00, ro = 27a 
and 77 ~ 0.7). The energies corresponding to bulk states are 
emphasized by purple shaded regions, see Fig. 2b. 
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evolution of this signal, which has non-vanishing values 
in the vicinity of the Fermi radius Rp, we show the chiral 
evolution of the edge states contribution p e d ge (x, t; ±$) 
in Figs. 8b-c. At small times, the overlap between 
the two contributions p e d ge (#, t\ ±$) decreases in time, 
leading to a progressive broadening of the signal 5p(x, t) 
along the ID circular edge. Then, after reaching a ro- 
tation of ~ 7r/4, the overlap increases, and eventu- 
ally leads to a vanishing of the signal 5p(x, t) « when 
the edge states have undergone a rotation of tt/2, where 
PedgeOM; +*) ~ Ped g e(x,t] -$). In Fig. 8, this happens 
at time t* ~ 49ft/ J, indicating that the edge states angu- 
lar velocity is ~ 0.03 J/h for Rp = 27a. The opposite- 
flux method therefore offers a general technique for em- 
phasizing the existence of chiral edge states in dispersive 
systems, and also, for evaluating their characteristic an- 
gular velocity. 

We verified that a slight difference in the filling, 
£f($+ = + 1 / 3 ) ~ #f($- = -1/3) ± 0.1J, or varia- 
tions in the flux, <£ + = 1/3 and ~ -$ + ± 0.01, does 
not significantly affect the signal 5p(x,t) shown in Fig. 
8 a, highlighting the robustness of this method against 
possible experimental imperfections. 
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FIG. 7: Wave-function overlaps between initial and fi- 
nal states, a, The overlap | (xa \<fi\) | 2 between the eigenstates 
of the initial (% a ) and final (0a) Hamiltonians, represented as 
a function of their energies. The parameters are the same as 
in Fig. 6. Note that the edge states are found in the bulk gap, 
namely, within the range E « [—2. 9 J, —1.3 J], b, Cut through 
the plot in Fig. a, for a specific initial edge state with energy 
E e = — 2J. The initial edge modes effectively project to the 
final edge modes, highly limiting the population of high en- 
ergy (dispersive) bulk states. The energies corresponding to 
bulk states are emphasized by purple shaded regions in a and 
b. 



FIG. 8: The opposite-flux method for dispersive 
systems, a, Evolution of the difference 5p = p(x,t;<f> = 
+ 1/3) — p(x,t;<& = —1/3), for the same configuration as in 
Figure 3. b, Evolution of the initially populated edge states 
pedge(x,t) for $ = +1/3, and c, for $ = —1/3. The vanishing 
of the signal 5p(x) ~ corresponds to a rotation of ir/2. This 
happens at time t* ~ 49ft/ J, indicating that the edge state 
angular velocity is ~ 0.03 J/h for Rf = 27a, in agreement 
with the results reported in Refs. [17]. 



measurements with opposite magnetic fluxes, is related 
to the edge states contributions p e dge(#, t; ±$). We show 
in Fig. 8a, the evolution of 5p(x, t) for the dispersive case 
$ = 1/3, which clearly indicates that the contributions 
from the bulk Pbuik(#,£; +$) ~ Pbuik(#,£; — vanish 
from the signal at all times t = 7 — 70ft/ J. Accordingly, 
Sp(x,t) w Ped g eO,£; +$)-ped g e(x,t; -$). To clarify the 



Appendix C: The edge-filter method for dispersive 
bulk bands 

Another strategy consists in allowing the edge states 
to propagate at t > 0, while forbidding the bulk states to 
penetrate the regions initially occupied by the holes. This 
can be achieved by suddenly lowering the additional po- 
tential walls Vhoie to some intermediate value V^oie ^ 0> 
instead of removing them completely at t = 0. This 
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"edge-filter" scheme is illustrated in Fig. 9a, for the dis- 
persive case $ = 1/3. By suddenly lowering the walls 
potential to the value ~ W, where W is the width 

of the lowest bulk band, we limit the undesired filling of 
the holes by the bulk states at times t > 0. In contrast, 
the populated edge states with energy e e > V^j^ + e m i n , 
where e m i n is the minimum of the bulk band, are allowed 
to propagate around the holes without being spoiled by 
the bulk. The resulting time evolution of the density 
p(x,t), presented in Fig. 9b, shows a clear propaga- 
tion of the edge states around the holes. Experimentally, 
this method offers an efficient method to isolate the edge 
states contribution from the spoiling bulk background, 
but it necessitates a very precise control over the poten- 
tial strength Vhoie- 
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FIG. 9: The edge-filter method for dispersive sys- 
tems, a, Energy spectrum E(k) as a function of the quasi- 
momentum k for <E> = 1/3. Also shown are the typical Fermi 
energy and the final walls potential strength V^J° ~ W, 
where W is the width of the lowest bulk band, b, Evolution 
of the spatial density p(x,t) for $ = 1/3, Ef = —0.85 J, and 
infinitely abrupt confinement /walls. At t = 0, the strength 
of the hole potential is suddenly reduced to = 0.5 J. 

Note that most of the bulk states are excluded from the holes 
during the evolution, yielding a clear edge-state signal, to be 
compared with Fig. 3 a. 



Appendix D: Sensitivity to imperfect filling 



the dispersionless case $ = 1/5. We remind that in this 
configuration, the clear separation of the bulk and edge 
states contributions to the density relies on the fact that 
the bulk states are dispersionless (they are described by 
a quasi-flat band). Here, we show that this picture holds 
even when the second (dispersive) bulk band is dramati- 
cally filled, see Fig. 10. From this result, we find that the 
high-energy bulk states contribute in a non-chiral manner 
to the density evolution, see Fig. 10, and that their dis- 
persive motion is slow compared to the edge states prop- 
agation along the circular boundary. In particular, this 
shows that the chiral picture drawn by the density p(x, £), 
and which will be imaged in an experiment, can unam- 
biguously be attributed to the populated edge states. We 
conclude that the edge-state signal obtained for the in- 
teresting case $ = 1/5 remains clear and detectable, as 
long as sufficiently many edge states are initially pop- 
ulated. In particular, this indicates that our scheme is 
robust at finite temperature T > 0, as long as it remains 
small compared to the gap's width A, in order to insure 
a sufficiently large edge-state population. 



1 ll|ii|fti! i ;r^r"'"V" i ; }non-chiral and dispersive bulk states 
J chiral states on the edge 
■ } non-chiral and dispersionless bulk states 
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FIG. 10: Sensitivity to imperfect filling in the flat- 
band configuration, a, Energy spectrum E(k) as a func- 
tion of the quasi- momentum k for <E> = 1/5, indicating the 
"dramatic" Fermi level Ef — —J used in b, as well as the 
characteristics of the occupied states, b, Evolution of the 
spatial density p(cc, t) for Ef — — J and infinitely abrupt con- 
finement/walls. Note the non-chiral behavior of the dispersive 
bulk states, to be compared with Fig. 4 a. 



Our general scheme is based on the possibility of 
preparing a QH atomic state, which can be achieved by 
generating a magnetic flux <I> in the lattice and filling 
the lowest bulk band completely. In other words, one 
has to tune the total number of particles N pa>rt in such a 
way that the Fermi energy E^ lies within the lowest bulk 
gap. According to the topological nature of the lowest 
bulk band, one is then guaranteed that topological edge 
states are populated. In practice, the total number of 
particles (and the corresponding Fermi energy E-p) can 
be tuned with a great precision in cold- atom experiments. 
However, it is instructive to test the robustness of our 
method against inexact filling effects, in particular, for 



Appendix E: Gaussian walls and smooth circular 
confinements 

When the magnetic flux $ = 1/g, with q G Z, the 
lowest bulk band is characterized by the Chern number 
v = — 1. This result, which only relies on the topological 
order associated with the bulk, guarantees the existence 
of a single edge mode whose energy lies within the 
first bulk gap. This "topological robustness" of the 
edge states is maintained when using smooth confining 
potentials and when adding walls Viioie to the system, as 
long as these features do not destroy the bulk gaps. In 
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this Section, we show that the edge states lying within 
the first bulk gap survive for realistic confinements and 
gaussian walls. We then demonstrate that the density 
evolution p(x,t), previously shown for infinitely abrupt 
potentials (main text), are qualitatively equivalent for 
smooth potentials. We finally comment on the fact that 
it is highly desirable to design sharp confining potentials 
to improve the edge state signal in experiments. 

In this Section, we focus on the case = 1/5, but we 
stress that the discussion holds for configurations that 
lead to similar band structures [16, 17, 26]). In the ab- 
sence of walls Vhoie = 0? the edge states lying in the first 
bulk gap are radially localized, with a radius determined 
by their energy and the circular confinement. Writing 
the circular confinement as V con f(r) = ^(r/ro) 7 , we find 
that an edge state <p e with energy e e is characterized by 
a localization radius 



7 = 10 



7 = 4 



n\M 2 



Re 



ro 



*\/Vo, 



(El) 



where e m i n is the minimum of the bulk band. This re- 
sult is illustrated in Fig. 11, for ro = 21a, Vq = J 
and 7 = oo,10,4, where the wavefunctions amplitudes 
\cj)\(x,y)\ 2 are plotted as a function of the x coordinate, 
and their corresponding energies e\. For an infinitely 
abrupt trap, 7 = 00, all the edge states are located 
at the constant Fermi radius Rp = ro- Therefore, the 
edge states contribution to the density p e dge yields a 
clear circular signal, with localization length of the or- 
der of the lattice spacing a. In contrast, for finite 7, 
the populated edge states are localized on different radii 
R e G [Rf — Sr, leading to a broadening of the edge- 
state signal Pedge- For the situation illustrated in Fig. 11, 
this broadening is of the order Sr ~ 5a for 7 = 10 and 
Sr ~ 10a for 7 = 4. Let us stress another crucial aspect 
of these smoothly confined QH systems, which is the fact 
that the angular velocity 6 of the chiral edge states, as 
well as the number N e d ge of available edge states within 
a bulk gap, highly depend on the potential's smoothness 
7: the angular velocity is maximized for highly abrupt 
confinements (7 — >• 00), while the number N e d ge is larger 
for smooth potentials [16, 17, 26]. We numerically eval- 
uated the angular velocity of the edge states for = 1/5 
and ro = 21a, and we found e ~ 0.06 J/ h for 7 = 00, 
e « 0.02 J/h for 7 = 10 and e < 0MJ/H for 7 = 4 (and 
we note that the angular velocity 6 e ex 1/R e ). Scaling to 
a more realistic radius ro = 100a, we find that the edge 
states, which are populated below Ep — —1.5 J, undergo 
a rotation of 7r/5 after a time 

t{6 = tt/5) - 50ft/ J, 7 = 00 (r = 100a), 

t{6 = tt/5) - 150ft/ ^ 7 = 10 (r = 100a), 

t(0 = tt/5) - 300ft/ J, 7 = 4 (r = 100a), 

indicating that it is highly desirable to design a sharp 
circular confining trap 7 ^> 10, to clearly observe the 
edge states rotating motion during reasonable experi- 
mental times t ~ 10 — 100ft/ J. 

We now investigate the density evolution p(x,i) for 
smooth confining traps and initial gaussian walls Vhoie- 
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FIG. 11: Smooth confinements and edge states. The 

amplitudes \(j)\(x, y)\ 2 of the single-particle wavefunctions as 
a function of the x coordinate and their corresponding energy 
e\, and setting y at the center of the trap. Here, $ = 1/5, 
and the external potential is given by V r con f(r) = Vb(r/ro) 7 , 
with ro — 21a and Vq = J. a: 7 = 00, b: 7 = 10, and 
c: 7 = 4. The Fermi radius Rf is shown for Ef — —1.5 J. 
The populated edge states <j) e are localized on the radii R e = 
R e (e e )i see Eq. (El), leading to a broadening 5r of their 
contribution to the spatial density p e dge for finite 7. 
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FIG. 12: Smooth confinements and the density evo- 
lution. Evolution of the spatial density for $ = 1/5 and 
Ef = —1.5 J. The walls Vhoie are produced by asymmetric 
gaussian potentials, with standard deviations a y w V2a x . 
The external confining potential is Vconf(r) = Vb(r/ro) 7 , with 
ro — 21a and Vq — J . a: 7 = 00, b: 7 = 10, and c: 7 = 4. 
Note the broadening Sr of the edge-state signal, see Fig. 11, 
and also, the deceleration of the motion as the confinement 
becomes smoother 7 = 00 — >• 4. 



First of all, we note that the presence of gaussian walls 
does not destroy the edge states lying within the first 
bulk gap. In this bat geometry, the edge states are 
localized on the outer edge delimited by V con f(r), but 
also, on the smooth boundary delimited by the gaussian 
potentials. Therefore, when 7 ~ 00, the edge states 
behave as in Fig. 11a in the vicinity of the outer 
circular edge r w ro, whereas they behave similarly as 
in Fig. 11c in the vicinity of the gaussian walls. We 
point out that, in our scheme, it is the behavior of the 
edge states near the outer circular edge which plays an 
important role. Indeed, as shown in Fig. 12a, replacing 
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the infinitely abrupt walls T4oie by gaussian potentials 
does not qualitatively affect the evolution of the density 
p(x, t) presented in Fig. 4a. 

However, as can be anticipated from the discussion 
above, replacing the perfectly sharp potential V con t(r) by 
smoother confinements, 7 = 00 — » 4, has dramatic con- 
sequences on the dynamics. In Fig. 12b, which shows 
the evolution of the density for 7 = 10, we clearly ob- 
serve the broadening Sr ~ 5a of the edge-state con- 
tribution p ec }ge as they progressively encircle the 
holes. We also note the slower motion undergone by 
the edge states, which have a reduced angular velocity 
#(7 = 10) ~ #(7 = oo)/3, see above. An even more dra- 
matic situation is illustrated for the case 7 = 4 in Fig. 
12c. These results demonstrate the robustness of the 
edge state motion in the presence of smooth confining 
traps and gaussian walls, but they greatly emphasize the 
importance of designing sharp external confinements to 
improve the experimental detectability of the topological 
edge states. 



Appendix F: The holes geometry and size effects 

Our setup illustrated in Fig. 1 features two large re- 
pulsive potentials, which create the initial bat geometry. 
These holes are chosen to be created by infinitely abrupt 
walls Vhoie? which are delimited by the two ellipses 

(x±r /2) 2 + (y/V2) 2 < r^/4, 

where the coordinates (x,y) = at the center of the 
trap. This choice is motivated by the fact that these walls 
coincide (up to first order) with the external circular wall 
V con f(r) of radius 7*0, in the vicinity of the poles located 
at (±7*0) 1<e- Note that in the following of this discussion, 
we consider that V con f(r) oc (r/ro) 7 with 7 = 00. 

The geometry produced by these potentials is particu- 
larly suitable to emphasize the edge states motion at time 
t > 0, as they encircle the regions initially surrounded by 
the elliptical walls. However, in this exotic geometry, the 
large holes correspond to a massive distortion of the stan- 
dard circular system. It is thus interesting to study the 
validity of our method, as we progressively reduce the 
size of the holes. Here, we demonstrate that our scheme 
still shows a clear manifestation of the edge states, even 
in the limit where the walls Vhoie on ly represent a small 
perturbation of the system, in the vicinity of its circular 
boundary. To analyze this, we propose to calculate the 
time evolution of the density p(x,t) after removing the 
elliptical walls delimited by the more general equations 



{x±(b-l)r /b) 2 +(y/Vb) 2 <{r /bf 



(Fl) 



which can be made arbitrarily small (b ^> 1), while main- 
taining the smooth contact with V con f(r) at the poles. 
The results are shown in Fig. 13 for the flat-band config- 
uration previously studied in Fig. 4 (for 6 = 2), but with 
smaller initial holes b = 4, 8, 16. This picture emphasizes 
the fact that our scheme still allows to detect the chiral 



motion of the edge states, in the limit of small perturba- 
tive walls Vhoie- However, we stress that it is crucial to 
prohibit any broadening of the edge state signal in the 
perturbative regime b ^> 2, which necessarily requires 
the use of an extremely abrupt external potential V con { 
with 7 ~ 00 (Appendix E). Moreover, considering smaller 
holes also demands to further reduce the bulk dispersion, 
which can be achieved by considering the quasi-flat-band 
configuration $ = 1/5, see Fig. 13, or by exploiting the 
"opposite-flux" or "edge-filter" methods. 
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FIG. 13: Effects of the holes size. The spatial density 
p(x,t) for $ - 1/5, E F = -1.5J, r - 27a and 7 = 00. The 
holes are described by Eq. (Fl) and correspond to a, b = 4, 
b, b = 8, and c, b = 16. Note that the chiral motion of the 
edge states is still visible, even in the limit of tiny perturbative 

Walls Vhoie- 



Appendix G: Robustness of the edge states against 
disorder 



We finally investigate the robustness of the edge states 
motion in the presence of disorder. This perturbation, 
which plays a fundamental role in solid-state systems, 
can be engineered in optical-lattice setups, e.g., using 
speckle potentials. We study the effects of disorder by 
considering a random potential Kand, with energies uni- 
formly distributed within the range V ran( ±(x) G [—D,D]. 
The results are presented in Fig. 14, for the case $ = 1/5. 
We find that the chiral edge states signal remains robust 
for disorder strengths D < 1.5 J ~ A, where A is the size 
of the bulk gap. Interestingly, we can still distinguish 
a cyclonic cloud - a signature of the edge states chiral- 
ity - for values up to D ~ 3 J ~ 2 A. The immunity of 
the chiral edge states against disorder, a hallmark of the 
quantum Hall effect, could thus be demonstrated using 
our cold-atom setup. 
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FIG. 14: Effects of disorder. The spatial density p(x,t) 
for $ = 1/5, Ef = —1.5 J, ro — 27a and 7 = 00. The disorder 
strength is a, D = 0.5J, b, D = 1J and c, D = 1.5J. 



